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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  provides  empirical 
equations  for  estimating  several  parameters  of  unsteady  flow  resulting  from  the  combination  of 
steady  storm  surge  overflow  and  overtopping  of  irregular  waves  at  a  trapezoidal-shaped  earthen 
levee.  Equations  are  given  for  the  average  overtopping  discharge  and  the  cumulative  probability 
distribution  of  instantaneous  overtopping  discharge.  On  the  landward-side  slope,  empirical  equations 
can  be  used  to  estimate  the  mean  flow  depth,  mean  flow  velocity,  root-mean-square  wave  height, 
and  velocity  associated  with  the  overtopping  wave  front.  Worked  examples  illustrate  application  of 
the  empirical  equations. 

INTRODUCTION  AND  BACKGROUND:  Earthen  levees  are  used  extensively  in  the  United 
States  to  protect  populations  and  infrastructure  from  periodic  floods  and  high  water  due  to  storm 
surges.  Overtopping  of  levees  and  dikes  produces  fast-flowing,  turbulent  water  velocities  on  the 
landward-side  slope  that  can  damage  the  protective  grass  covering  and  expose  the  underlying  soil  to 
erosion.  If  overtopping  continues  long  enough,  the  erosion  may  eventually  result  in  loss  of  levee 
crest  elevation  and  perhaps  breaching  of  the  protective  structure. 

Economics  often  dictate  levee  designs  with  crown  elevations  having  a  risk  that  some  wave/surge 
overtopping  will  occur  during  extreme  events.  In  addition,  crown  elevations  for  older  levee  systems 
may  have  been  established  without  complete  information  about  possible  water  levels  that  might 
occur  during  extreme  events.  Even  levees  that  presently  have  sufficient  freeboard  to  withstand  all  but 
the  most  extreme  storms  may  become  vulnerable  to  wave  overtopping  and  storm  surge  overflow  in 
the  future  if  sea  level  continues  to  rise  at  projected  rates.  The  cost  and  engineering  challenges  of 
raising  all  levees  to  elevations  where  overtopping  will  be  within  tolerable  limits  may  be 
insurmountable.  Levees  that  cannot  be  raised  remain  vulnerable,  and  the  landward-side  levee  slopes 
must  be  protected  with  some  type  of  strengthening  alternative  such  as  turf  reinforcement,  soil 
strengthening,  or  hard  armoring. 

Assessment  of  levee  reliability  and  design  of  slope  protection  alternatives  requires  estimates  of  wave 
overtopping  and  stonn  surge  overflow  that  will  occur  for  a  specified  set  of  storm  parameters.  This 
CHETN  provides  overtopping  estimation  equations  for  the  case  of  steady  overflow  due  to  a  still 
water  level  above  the  levee  crown  elevation  combined  with  irregular  wave  overtopping.  The 
overtopping  flow  is  unsteady  in  time  and  spatially  nonuniform,  as  illustrated  in  Figure  1.  As  each 
overtopping  wave  passes  over  the  levee  crown,  it  has  a  somewhat  triangular-shaped  discharge 
distribution  with  a  maximum  discharge  at  the  leading  edge  that  is  several  times  greater  than  the  time- 
averaged  discharge.  If  the  storm  surge  elevation  above  the  levee  crown  is  small  and/or  the  waves  are 
large,  the  levee  crown  and  landward-side  slope  will  “go  dry”  during  wave  troughs.  This  fluctuating 
discharge  may  contribute  to  accelerated  erosion  of  the  landward-side  levee  slope. 
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Figure  1.  Sequence  of  wave  overtopping  on  a  scale-model  levee. 


LABORATORY  EXPERIMENTS:  Laboratory  experiments  of  wave  overtopping  combined  with 
steady  surge  overflow  were  conducted  at  a  nominal  prototype-to-model  length  scale  of  25-to-l  in  a 
3-ft-wide  wave  flume  at  the  U.S.  Army  Engineer  Research  and  Development  Center’s  Coastal  and 
Hydraulics  Laboratory.  The  tested  levee  cross  section  replicated  in  the  physical  model  is  shown  in 
Figure  2.  Seaward  is  on  the  left  side  of  the  figure,  and  the  dimensions  correspond  to  full-scale  units. 
This  cross  section  was  typical  of  the  Mississippi  River  Gulf  Outlet  that  experienced  severe 
overtopping  during  Hurricane  Katrina. 
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Figure  2.  Levee  cross  section  tested  in  physical  model  (full-scale  dimensions). 


Prototype-scale  target  wave  and  surge  parameters  for  testing  were  three  surge  elevations  (5  =  + 1 .0, 
+3.0,  and  +5.0  ft  above  levee  crest),  three  significant  wave  heights  (Hm o  =  3.0,  6.0,  and  9.0  ft),  and 
three  peak  wave  periods  (Tp  =  6,  10,  and  14  sec).  This  gave  a  total  of  27  unique  conditions  for 
combined  wave  and  surge  overtopping. 

Instantaneous  discharge  over  the  levee  crest  was  estimated  at  a  location  near  the  landward-side  edge 
of  the  levee  crest.  Coincident  values  of  water  depth  (measured  with  a  pressure  gauge)  and  horizontal 
velocity  (measured  with  a  laser  Doppler  velocimeter)  were  multiplied  at  each  instant  in  time  to 
determine  the  time  series  of  instantaneous  discharge.  This  estimate  assumed  that  velocity  was 
horizontal  and  constant  throughout  the  water  column  at  that  location.  Incident  wave  characteristics 
were  measured  at  a  three-gauge  array  located  seaward  of  the  levee. 

For  overtopping,  the  freeboard  is  defined  as  the  difference  between  the  levee  crest  elevation,  hc,  and 
the  still- water  elevation,  hs,  i.e.,  Rc  =  (hc  -  hs).  When  the  water  elevation  is  higher  than  the  levee 
crest,  Rc  <  0,  which  is  referred  to  as  negative  freeboard. 

AVERAGE  DISCHARGE  FOR  COMBINED  WAVE  AND  SURGE  OVERTOPPING:  For 

each  of  the  27  experiments,  an  average  was  calculated  over  the  portion  of  the  instantaneous 
discharge  time  series  that  included  both  wave  overtopping  and  steady  overflow.  This  average 
represents  the  average  overtopping  discharge  for  that  particular  combination  of  negative  freeboard 
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and  incident  irregular  waves.  It  is  customary  in  overtopping  research  to  express  the  average  wave 
overtopping  discharge  in  terms  of  a  dimensionless  parameter  defined  as 


Qm 


Vs  Hl« 


(i) 


where 


Qws  =  dimensionless  average  discharge  due  to  combined  waves  and  surge 
qws  =  average  discharge  due  to  combined  waves  and  surge  per  unit  crest  length 
g  =  acceleration  of  gravity 
Hm o  =  energy-based  significant  wave  height 

Figure  3  plots  the  dimensionless  combined  wave/surge  average  overtopping  discharge  versus  the 
relative  (negative)  freeboard  for  all  27  experiments.  The  indicated  surge  levels  in  the  plot  legend  are 
the  average  of  the  negative  freeboards  determined  for  all  nine  experiments  at  each  nominal  surge 
level. 
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Figure  3.  Dimensionless  combined  wave/surge  average  discharge  versus 

relative  freeboard. 


As  shown  in  Figure  3,  the  measurements  gave  a  nice  trend  with  increasing  relative  freeboard,  and  the 
solid  line  is  a  best-fit  empirical  equation  given  by  the  formula 


3 


ERDC/CHL  CHETN-lll-78 
May  2008 


Qws 


?w, 

4  s  Hl  0 


0.0336  +  0.53 


/  \  1.58 

-R  ' 


V  y 


(2) 


This  best- fit  equation  had  a  correlation  coefficient  of 0.9987  and  a  root-mean-square  (RMS)  percent 
error  of  0.12.  (Note  that  /?,  must  be  specified  as  a  negative  number  so  the  ratio  in  brackets  will  be 
positive.)  Peak  spectral  wave  period  had  negligible  influence  in  the  determination  of  Qws  for  the 
range  of  periods  tested  in  the  model. 

Like  any  empirical  equation,  application  of  Equation  2  should  be  limited  to  the  range  of  tested 
parameters.  In  particular,  seaward  (flood-side)  levee  slopes  different  from  1 :4.25  could  influence  the 
wave  overtopping,  but  seaward  slope  effects  should  decrease  as  surge  level  increases.  For  zero 
freeboard  ( Rc  =  0),  Equation  2  yields  a  constant  dimensionless  average  wave  overtopping  rate  that 
falls  near  the  middle  of  the  range  of  experiment  results  given  by  Schiittrumpf  et  al.  (2001).  This  is 
encouraging;  however,  estimates  for  zero  freeboard  should  be  made  using  Schiittrumpf  s  equations 
as  presented  in  the  Overtopping  Manual  (Pullen  et  al.  2007)  rather  than  Equation  2. 

PROBABILITY  DISTRIBUTION  OF  INSTANTANEOUS  DISCHARGE:  During  combined 
wave  and  surge  overtopping,  each  wave  will  have  instantaneous  peak  discharge  that  can  be  several 
times  the  average  overtopping  rate,  qws.  These  peak  discharges,  while  short  in  duration,  may  be 
responsible  for  much  of  the  erosion  of  levee  soil  or  instability/failure  of  annoring  alternatives.  Thus, 
prediction  of  the  instantaneous  overtopping  distribution  is  needed  for  more  effective  design  of  levee 
slope  protection  alternatives. 

The  time  series  measurements  of  overtopping  flow  discharge  acquired  near  the  landward-side  edge 
of  the  levee  crest  were  analyzed  to  determine  the  cumulative  distribution  of  instantaneous  discharge 
over  the  levee  crest.  The  analysis  revealed  that  still-water  elevation  above  the  levee  crest  (i.e.,  the 
steady  overflow  discharge)  is  the  most  important  hydrodynamic  parameter,  and  it  most  closely 
controls  the  scale  of  the  cumulative  probability  distribution.  Wave  height  and  period  appear  to 
influence  the  shape  of  the  distribution  extreme  tail. 

The  measured  cumulative  discharge  distributions  were  fit  to  the  Weibull  cumulative  probability 
distribution  given  by  the  expression 


P(q  <  q„)  =  1-exp 


'  qA 

b 

l  c  J 

(3) 


where  c  is  the  “scale  factor”  and  b  is  the  “shape  factor”  of  the  distribution.  The  scale  factor,  c,  has 
units  of  discharge  per  unit  length,  whereas  the  shape  factor,  b,  is  dimensionless.  Equation  3  gives  the 
probability  that  a  value  of  discharge  will  be  below  the  specified  discharge,  q*.  The  corresponding 
distribution  of  percent  exceedance  is  given  by 


100 ■ exp 

r  qAb 

l  C  J 

(4) 
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The  best-fit  procedure  produced  values  of  the  distribution  parameters,  b  and  c.  Figure  4  shows  an 
example  of  a  good  fit  (upper  plot)  and  a  mediocre  fit  (lower  plot)  of  the  Weibull  cumulative 
distribution  to  the  measured  distributions  (given  as  percent  exceedance  distributions).  The  curves  are 
plotted  with  a  semi-log  ordinate  to  better  present  the  extreme  end  of  the  exceedance  probability. 
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Figure  4.  Example  best  fits  of  Weibull  cumulative  probability  distribution  to  instantaneous  discharge. 


The  curves  in  Figure  4  give  the  percent  occurrence  that  the  overtopping  discharge  will  be  above  a 
given  value.  For  example,  in  the  upper  plot  of  Figure  4,  for  10  percent  of  the  time  the  discharge  will 
be  above  q~  23  ft3/sec  per  foot,  and  the  discharge  will  exceed  q~  36  ft3/sec  per  foot  for  1  percent  of 
the  time.  The  average  wave/surge  overtopping  for  this  experiment  was  qws=  13.1  ft3/sec  per  foot. 
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The  dimensionless  shape  factor,  b,  exerts  control  over  the  extreme  tail  of  the  distribution,  and  the  tail 
is  sensitive  to  small  differences  in  the  shape  factor.  Through  inspection,  and  after  several  attempts 
using  dimensionless  combinations  of  the  wave  height,  wave  period,  and  steady  surge  discharge,  the 
best  empirical  result  for  the  shape  factor  was  the  relationship 


b  =  8.10 


s  X  0.343 

<ls 

ygH,nOTP , 


(5) 


where  qs  is  estimated  from  the  broad-crest  weir  formulation  as 


<ls  = 


0.5443^  |^|3'2 


(6) 


Figure  5  shows  the  best-fit  curve  for  the  shape  factor  given  by  Equation  5  along  with  the  plotted 
values  of  b  obtained  from  the  measured  data.  The  correlation  coefficient  for  the  best  fit  was  0.922. 


The  three  data  points  in  Figure  5  with  values  of  b  greater  than  3  deviated  significantly  from  the  data 
trend.  These  points  came  from  experiments  with  the  surge  level  at  approximately  +4.3  ft  and 
nominal  wave  heights  around  3  ft  (scaled  from  model  to  prototype).  A  possible  explanation  for  the 
observed  deviation  might  be  the  effect  of  wave  deformation  by  the  strong  steady  overflow  as  waves 
reach  the  levee.  The  nonlinear  interaction  between  the  waves  and  the  coincident  current  at  the  levee 
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crest  is  more  pronounced  at  the  higher  surge  level,  and  smaller  waves  would  be  affected  more  than 
larger  waves.  The  resulting  smaller  wave  heights  would  reduce  the  extreme  values  of  overtopping 
discharge,  giving  higher  values  for  the  shape  factor,  b. 


The  mean  of  the  Weibull  cumulative  distribution  is  equal  to  the  average  overtopping  discharge,  qws, 
and  it  can  also  be  expressed  in  tenns  of  the  b  and  c  parameters  as  (e.g.,  Goda  2000) 


<G, 


=  c  ■ 


1  + 


n 

b; 


(7) 


where  T  is  the  mathematical  gamma  function.  Equation  7  was  used  to  estimate  qws  using  the  best- fit 
values  of  the  Weibull  b  and  c  parameters  for  all  experiments,  and  comparison  with  values  of  qws 
determined  directly  from  the  discharge  time  series  showed  excellent  agreement.  Therefore, 
Equation  7  can  be  rearranged  to  give  an  expression  for  the  scale  factor,  c,  i.e., 


c  = 


1  + 


r 

b) 


(8) 


The  gamma  function  can  be  accurately  calculated  over  the  range  of  shape  factor,  b,  found  for  these 
experiments  (0.5  <  b  <  4)  by  the  best  fit  of  a  cubic  equation  given  by 

f(x)  =  1 .3  -  0.368  ■  x  -  0.072  ■  x2  +  0.09  ■  x3  (9) 


which  is  valid  for  1.25  <  x  <  3.  This  approximation  can  be  used  for  calculations  when  the  gamma 
function  is  not  available. 


In  summary,  the  cumulative  probability  distribution  of  instantaneous  discharge  over  the  levee  crest 
due  to  combined  wave  and  surge  overtopping  can  be  estimated  using  Equation  3  or  4  with  the 
Weibull  shape  factor,  b,  determined  from  Equation  5,  and  the  Weibull  scale  factor,  c,  estimated 
using  Equation  8  as  a  function  of  qws  (from  Equation  2).  For  an  application  in  which  calculation  of 
the  gamma  function  is  not  available,  use  Equation  9  within  the  stated  range  of  applicability. 

FLOW  PARAMETERS  ON  THE  LANDWARD-SIDE  LEVEE  SLOPE:  Steady  overflow  on  the 
steep  landward-side  slope  of  a  levee  is  supercritical  with  slope-parallel  velocities  increasing  down 
the  slope  until  a  balance  is  reached  between  the  momentum  of  the  flow  and  the  frictional  resistance 
force  of  the  slope  surface.  Once  tenninal  velocity  is  reached  (mean  velocity  is  constant),  the  steady 
flow  can  be  analyzed  using  the  Chezy  or  Manning  equation.  Flow  down  the  landward-side  slope 
caused  by  combined  waves  and  surge  overtopping  is  unsteady  and  more  difficult  to  analyze 
thoroughly.  In  this  section  empirical  equations  are  given  to  characterize  several  representative 
parameters  of  wave-related  unsteady  flow  on  the  landward-side  levee  slope. 

Average  Flow  Depth  on  Landward-Side  Slope.  The  average  flow  depth,  rjm,  perpendicular  to 
the  slope  was  calculated  at  four  pressure  gauge  locations  that  were  equally  spaced  down  the 
landward-side  slope.  Comparison  among  the  four  gauges  indicated  the  values  of  rjm  were  nearly 
constant  for  each  experiment,  so  an  average  was  taken  to  represent  the  mean  flow  thickness  on  the 
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landward-side  slope.  A  correlation  was  sought  between  rjm  and  the  hydrodynamic  forcing,  and  a 
reasonable  result  is  shown  in  Figure  6.  Values  have  been  scaled  to  prototype  dimensions  using  a 
length  scale  ratio  of  25-to- 1 .  The  solid  line  is  a  linear  best- fit  equation  given  by  the  simple  empirical 
expression 


yfgtf, 


2.226 


(10) 


This  equation  had  a  correlation  coefficient  of  0.988,  and  an  RMS  percent  error  of  0.134. 


Figure  6.  Average  overtopping  discharge  versus  mean  flow  thickness  on  landward-side 

slope  (prototype  scale). 


If  the  mean  velocity  on  the  landward-side  slope  is  defined  as  vm  =  qws  /  r/m,  then 

v,„=  2.226  (11) 

which  bears  resemblance  to  the  Chezy  equation. 

The  Chezy  equation  for  wide  channels  (hydraulic  radius  is  approximately  equal  to  flow  depth,  d  ) 
and  steady  flow  (friction  slope  is  equal  to  landward-side  levee  slope,  9)  is  given  as 
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v  =  ]^p'V^  (12) 

for  steep  slopes.  In  Equation  12  the  Chezy  coefficient  was  replaced  with  a  function  of  the  Fanning 
friction  factor  given  by  the  expression 


Under  the  assumption  the  Chezy  equation  is  an  appropriate  model  to  approximate  the  average  of  this 
rapidly  varying  unsteady  flow  situation,  it  is  hypothesized  that  the  constant  in  Equations  1 0  and  1 1  is 
a  function  of  both  levee  slope  and  a  representative  friction  factor,  as  given  by  the  first  radical  tenn  in 
Equation  12. 

The  Chezy  coefficient  Cz  varies  with  flow  thickness  on  the  slope;  therefore,  the  friction  factor,//-, 
also  varies  continuously  during  the  unsteady  overtopping  flow.  A  representative  friction  factor  for 
these  experiments  was  found  by  assuming  that  the  constant  in  Equation  1 1  is  equal  to  the  first  radical 
tenn  in  Equation  12.  This  yielded  a  value  offF=  0. 128.  However,  at  this  point  there  is  no  method  for 
estimating  similar  representative  friction  factors  for  other  slope  roughness,  so  little  benefit  is  gained 
by  including  friction  factor  at  this  time.  It  is  assumed,  however,  that  the  representation  of  structure 
slope  in  Equation  12  is  reasonable.  For  the  l-on-3  slope  in  these  experiments,  the  constant  in 
Equations  10  and  1 1  was  detennined  to  be  2.226  =  3.96  ^sinO  . 

Substituting  for  the  constant  in  Equation  10  and  rearranging  gives  a  tentative  equation  for  the  mean 
flow  depth,  i.e., 


nm  =  0-4 


-|l/3 


gsin# 


(tfvJ 


2/3 


and  the  mean  flow  velocity  equation  becomes  the  following 

v,n  = 2-5  {qm-g-  sin0)1/3 


(14) 


(15) 


The  constants  in  Equations  14  and  15  are  related  to  the  slope  roughness,  and  the  equations  are 
applicable  only  for  landward-side  slopes  having  roughness  similar  to  that  of  the  laboratory  experi¬ 
ments.  Friction  factors  for  grass  levee  slopes  should  not  be  too  much  higher  than  those  of  the 
experiments,  but  slopes  armored  with  riprap  or  similar  material  will  have  significantly  higher 
representative  friction  factors.  More  work  is  needed  to  determine  appropriate  representative  friction 
factors  for  a  range  of  slope  roughness. 

Estimation  of  Hrms  on  the  Landward-Side  Slope.  Values  of  root-mean-square  wave  height, 
Hrms,  perpendicular  to  the  slope  at  the  four  pressure  gauge  locations  on  the  landward-side  slope  were 
averaged  based  on  the  observation  that  little  variation  was  shown  over  this  down-slope  distance.  A 
relationship  was  sought  that  expressed  Hrms  in  tenns  of  other  parameters  that  could  be  estimated. 
Figure  7  shows  the  best  correlation  of  those  attempted.  Wave  period  was  found  to  have  only 
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marginal  influence,  and  some  of  the  observed  scatter  is  possibly  explained  by  wave  period  variation. 
The  solid  curve  represents  the  best  fit  of  a  one-parameter  exponential  function  given  as 


H„ 


Vn 


-  3.43 ■ exp 


V  ^ m0  J 


(16) 


This  best-fit  equation  had  a  correlation  coefficient  of  0.966  and  an  RMS  percent  error  of  0.094. 
When  applying  Equation  16,  keep  in  mind  that  Rc  must  be  entered  as  a  negative  number,  and  the 
equation  should  not  be  applied  for  cases  in  which  Rc  >  0  (surge  level  is  at  or  below  the  crest 
elevation). 


R  /  H 

c 


mO 


Figure  7.  Estimation  of  Hrms  on  the  landward-side  slope  as  a  function  of  r/m. 


The  effect  of  levee  landward-side  slope  is  included  in  Equation  16  through  the  estimation  of  the 
mean  flow  depth,  using  Equation  14.  However,  this  is  still  tentative  and  may  not  be  correct,  so 
caution  is  urged  until  additional  verification  of  Equation  16  is  forthcoming.  Also,  note  that  HrmslTjm 
is  less  than  unity  when  Rc/Hmo  is  less  than  -1.2.  In  these  cases,  the  landward-side  slope  does  not  go 
dry  during  passage  of  the  wave  trough,  and  the  mean  flow  depth  is  greater  than  the  crest-to-trough 
length  represented  by  Hrms. 
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Once  Hrms  is  estimated,  the  Rayleigh  distribution  can  be  used  to  obtain  reasonable  estimates  for 
and  //mo  according  to  the  equations 

HU3=IM6-Hrms  and  Hmo  =\.80- Hrms  (17) 

Comparisons  between  measured  characteristic  wave  heights  and  those  predicted  by  the  Rayleigh 
distribution  using  measured  values  of  Hrms  were  performed  for  all  laboratory  data.  An  excellent 
prediction  was  found  for  Hy^,  and  the  prediction  lor  //mo  was  quite  good  with  a  little  more  scatter 
and  a  tendency  to  overpredict  slightly.  However,  the  Rayleigh  distribution  overpredicted  the 
measured  values  of  T/i/ioo  by  as  much  as  25  to  40  percent,  indicating  that  the  highest  waves  do  not 
conform  well  to  the  Rayleigh  distribution.  Nevertheless,  the  reasonable  fit  of  the  Rayleigh 
distribution  is  a  useful  finding  because  it  allows  much  of  the  wave  height  distribution  on  the 
landward-side  slope  to  be  characterized  in  tenns  of  the  RMS  wave  height,  Hrms. 

Estimation  of  Wave  Front  Velocity  on  the  Landward-Side  Slope.  Flow  velocity  was  not 
measured  directly  on  the  landward-side  slope  during  the  laboratory  experiments.  However,  it  was 
possible  to  make  a  rough  estimate  of  the  speed  of  the  wave  front  by  dividing  the  distance  between 
the  first  and  last  pressure  gauges  by  the  time  it  took  the  wave  front  to  move  over  that  distance.  This 
approximate  wave  front  velocity  represents  an  average  over  the  distance  between  the  four  pressure 
gauges.  In  fact,  the  flow  might  have  been  accelerating  over  this  reach.  The  accuracy  of  wave  front 
velocity  estimates  using  this  technique  is  limited  by  the  data  sampling  rate  of  50  Hz  used  in  the 
experiments.  One  data  point  shift  either  way  is  a  relatively  large  change  in  estimated  velocity.  For 
this  reason,  the  wave  front  velocity  analysis  was  restricted  to  the  first  and  last  gauges  to  maximize 
the  distance  between  gauges,  and  thus,  to  reduce  the  potential  for  error. 

Figure  8  plots  the  wave  front  velocity  versus  the  velocity  of  a  bore  having  a  frontal  depth  equal  to 
Hrms.  Values  of  Hrms  were  the  averages  taken  over  the  four  pressure  gauges.  The  discrete  jumps  in 
velocity  on  the  plot  ordinate  represent  time  shifts  of  four,  five,  six,  and  seven  time  steps  (highest  to 
lowest  velocity).  The  actual  wave  front  velocities  of  some  points  are  probably  somewhere  between 
the  levels,  and  this  could  reduce  the  scatter  somewhat.  The  straight  line  shown  in  Figure  8  is  the  best 
fit  given  by  the  equation 


v„=3.85  (18) 

where  vw  is  the  wave  front  velocity.  This  best  fit  had  a  correlation  coefficient  of  0.89 1  and  an  RMS 
percent  error  of  0.114.  The  coefficient  in  Equation  18  is  constant  for  this  particular  data  set; 
however,  it  may  be  a  function  of  slope  angle  and  surface  roughness  even  though  Equation  16  for 
estimating  Hnns  explicitly  includes  slope  and  implicitly  includes  friction  effects  for  smooth  slopes. 


11 


ERDC/CHL  CHETN-lll-78 
May  2008 


Figure  8.  Estimated  wave  front  velocity  determined  from  the  first  and  last  pressures  gauges. 


Example:  Combined  Wave/Surge  Overtopping  Flow  Parameters 


Find:  The  flow  parameters  associated  with  wave  overtopping  combined  with  overflow  by  a  surge 
elevation  that  is  1 .3  ft  above  the  levee  crest  elevation.  The  levee  surface  is  considered  smooth  with 
roughness  comparable  to  grass  or  turf  reinforcement  mats. 

Given: 


Hm  o 

TP 

hc 

hs 

Rc 
tan  a 
tan  0 
g 


6  ft 

=  zeroth-moment  significant  wave  height 

8  sec 

=  wave  period  associated  with  the  spectral  peak 

17.0  ft 

=  levee  crest  elevation 

18.3  ft 

=  storm  surge  elevation 

-1.3  ft 

=  freeboard  [=  hc  -  hs\ 

1/4 

=  seaward-side  levee  slope 

1/3 

=  landward-side  levee  slope  [0  = 

18.4  deg] 

32.2  ft/sec2 

=  acceleration  due  to  gravity 
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The  seaward-side  levee  slope  is  not  used  in  these  calculations,  but  it  is  wise  to  make  sure  the  slope  is 
not  too  different  from  the  l-on-4.25  seaward-side  slope  used  in  the  physical  model  study  on  which 
these  equations  are  based. 

Calculate  Average  Overtopping  Discharge.  The  dimensionless  average  overtopping  discharge 
is  found  using  Equation  2. 


0, 


0.0336  +  0.53 

0.0336  +  0.53 


{Hm  0) 

-  (—1.3  ft) 
,  6ft 


\  1.58 

J 


a, =o.o8i 


And,  the  dimensional  overtopping  is  calculated  by  rearranging  Equation  1  to  yield 

=  a/777  Q„  =  y(32.2ft/secJ)(6ft’)  0.081 
qws  =6.75ft2/sec 

(Note  the  discharge  units  square  feet  per  sec  are  equivalent  to  cubic  feet  per  sec  per  foot,  or  cfs/ft.) 

Calculate  Instantaneous  Overtopping  Discharge  Cumulative  Distribution.  The  Weibull 
percent  exceedance  cumulative  probability  distribution  is  given  by  Equation  4.  The  first  step  in 
determining  the  distribution  parameters  is  estimating  the  steady  overflow  discharge  (Equation  6) 
due  to  the  surge  elevation  exceeding  the  crest  elevation. 

q,  =0.5443-^ |f?c|372  =  0.5443^32.2 ft/sec2  |-1.3ft|3/2 
qs  =  4.58ft2/sec 


(Note  that  qws/qs  =  1.47.)  Now,  calculate  the  Weibull  distribution  shape  factor,  b,  using  Equation  5. 

4.58ft2/sec 


,  0.343 


6  =  8.10 


7v 

y  &  ^m0  ^ p  J 


,  0.343 


=  8.10 


(32.2  ft/sec2)  (6  ft)  (8  sec) 


=  1.1 


The  Weibull  scale  factor,  c,  is  given  by  Equation  8.  For  convenience,  evaluate  the  gamma  function 
in  Equation  8  using  Equation  9  withx  =  1  +  Mb  =  1.91,  i.e., 

T(x)  =  1.3-  0.368  ■  x  -  0.072  ■  x2  +  0.09  ■  x3 

T(1 .9 1)  =  1 .3  -  0.368  ■  (1 .9 1)  -  0.072  ■  (1 .9 1)2  +  0.09  ■  (1 .9 1)3 

r(1.91)  =  0.962 
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Substituting  for  the  gamma  function  in  Equation  8  yields 

qws  _6.75ft2/sec  6.75ft2/sec 


r  n 

r  i  'i 

0.962 

f 

1  +  - 

f 

i+ — 

V  b) 

V  1.1/ 

7.02  ft  /sec 


The  percent  exceedance  cumulative  probability  distribution  for  this  example  becomes 

>1-1' 


P%  ( q  >  q ,)  =  100  ■  exp 


q * 


7.02  ft  /sec 


As  an  example,  the  percentage  of  instantaneous  overtopping  discharge  exceeding  q*  =  15  cfc/ft  is 
given  by 


P%  (q  >  15ft2/sec)  =  100  exp 


15ft2/sec 

v7.02ft2/secy 


10.0% 


In  other  words,  10  percent  of  the  instantaneous  discharge  values  are  greater  than  2.2  times  the 
average  overtopping  discharge. 

Calculate  Flow  Parameters  on  the  Landward-Side  Levee  Slope.  Under  the  assumption 
that  the  levee  surface  roughness  is  small  (similar  to  the  smooth  slope  used  in  the  laboratory  tests), 
the  mean  flow  thickness  perpendicular  to  the  landward-side  slope  is  estimated  from  Equation  14  as 


nm  =0.4 


-|l/3 


gsint? 


(a»J 


,2/3 


■  0.4 


1 


-|l/3 


(32.2ft/sec2)sin(18.4°) 


(6.75ft2/sec)2/3 


nm  =  0.4  (o.462s2/3/ft1/3  )(3 .572  ft4/3/s2/3 )  =  0.66  ft 


The  mean  flow  velocity  on  the  landward-side  slope  can  be  estimated  using  Equation  1 5 ,  but  an  easier 
method  is  using  the  definition 


v_  = 


_  l  ws  _ 


6.75ft2/sec 


=  10.2  ft/sec 


nm  0.66  ft 

The  RMS  wave  height  (perpendicular  to  the  levee  slope)  is  estimated  from  Equation  16  as 


H, 

H 


rms  i  m 


3.43 ■ exp 


r  R  A 


V  Hm«  J 


=  (0.66  ft) 


3.43  ■  exp 


^-1.3  fO 
6  ft 


1.8  ft 
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Because  the  wave  heights  on  the  landward-side  slope  appear  to  follow  the  Rayleigh  distribution 
except  at  the  extreme  tail,  other  representative  wave  heights  can  be  estimated  from  Equation  17,  i.e., 

HU3  =  1.416  •  Hrms  =  1.416  •  (1.8ft)  =  2.6ft 
Hino  =  1-80  ■  Hrms  =  1.80  ■  (1.8ft)  =  3.2  ft 


Finally,  a  rough  estimate  of  the  velocity  of  the  overtopping  wave  front  on  the  landward-side  slope  is 
given  by  Equation  1 8  as 


vlv  =3.85  ^gllrms  =  3.85  V(32.2ft/sec2)(1.8ft) 
vw  =  29.3ft/sec 


The  velocity  of  the  wave  front  is  almost  3  times  greater  than  the  mean  velocity  on  the  landward-side 
slope.  Most  importantly,  the  wave  front  velocity  estimate  should  be  considered  preliminary,  and  the 
margin  for  error  is  significant  due  to  the  method  by  which  the  equation  was  developed. 

Flow  Parameter  Variation  for  Same  Average  Overtopping  Discharge.  Different  combi¬ 
nations  of  incident  wave  height  and  freeboard  can  produce  the  same  value  of  average  overtopping 
discharge,  qws.  Table  1  shows  overtopping  flow  parameters  calculated  for  three  additional  cases  that 
produced  the  same  average  discharge  (qws  =  6.75  ft2/sec)  detennined  above. 


Table  1 

Additional  Parameter  Estimates  for  the  Same  Value  of  qws 

HmO 

Rc 

Qs 

c 

Tjm 

Vm 

Hrms 

vw 

(ft) 

(ft) 

ISBM 

(ft2/sec) 

b 

(ft2/sec) 

(ft) 

(ft/sec) 

(ft) 

(ft/sec) 

1.0 

-1.64 

6.75 

6.49 

2.3 

7.60 

0.66 

10.2 

0.44 

14.5 

3.0 

-1.59 

6.75 

6.22 

1.5 

7.50 

0.66 

10.2 

1.33 

25.2 

6.0 

-1.30 

6.75 

4.58 

1.1 

7.02 

0.66 

10.2 

1.80 

29.3 

9.0 

-0.75 

6.75 

2.00 

0.7 

5.46 

0.66 

10.2 

2.08 

31.5 

REMARKS:  All  four  of  the  cases  given  in  the  table  above  have  different  combinations  of  incident 
significant  wave  height  and  levee  freeboard,  although  the  average  overtopping  discharge,  qws,  was 
the  same.  The  mean  flow  parameters  on  the  landward-side  slope  ( ijm  and  vm)  were  also  identical. 
However,  the  parameters  of  the  Weibull  distribution  of  instantaneous  discharge  ( b  and  c),  the  RMS 
wave  height  (Ht7ns),  and  the  velocity  of  the  wave  front  (vw)  were  quite  different  among  the  cases.  The 
difference  of  wave  height  and  wave  front  velocity  on  the  landward-side  slope  suggests  that 
characterizing  levee  safety  in  terms  of  an  upper  limit  for  average  overtopping  discharge  may  not  be 
prudent  unless  it  can  be  shown  that  levee  erosion  and  armor  stability  is  more  a  function  of  the  mean 
flow  characteristics  than  a  function  of  the  variations  in  the  unsteady  overtopping  flow. 

SUMMARY:  This  CHETN  summarized  new  equations  for  average  overtopping  discharge  and 
distribution  of  instantaneous  discharge  associated  with  combined  wave  and  storm  surge  overtopping 
of  levees.  Equations  are  also  given  for  mean  flow  thickness,  RMS  wave  height,  mean  velocity,  and 
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velocity  of  the  wave  front  down  the  landward-side  slope.  An  example  problem  illustrates  application 
of  the  new  equations.  These  equations  are  intended  for  use  in  preliminary  design. 

The  empirical  equations  are  based  on  small-scale  laboratory  experiments  featuring  a  levee  with  a 
seaward-side  slope  of  1:4.25,  and  the  equations  may  not  be  applicable  for  levees  having  different 
seaward-side  slopes.  It  is  hypothesized  that  seaward-side  slope  may  not  be  as  important  for 
combined  wave  and  surge  overtopping  as  it  is  for  wave-only  overtopping.  Additionally,  several 
assumptions  were  invoked  during  equation  development.  First,  it  was  assumed  that  unsteady 
overtopping  flow  velocities  over  the  levee  crest  are  horizontal  and  unifonn  throughout  the  depth  at 
each  instant  in  time.  Second,  the  equations  for  the  landward-side  slope  include  slope  angle  and  a 
constant  that  is  likely  a  function  of  some  representative  friction  factor  characterizing  slope 
roughness.  Applicability  of  these  equations  for  landward-side  slopes  different  from  l-on-3  is 
uncertain,  and  the  equations  will  give  incorrect  estimates  where  slope  surface  roughness  is  not 
relatively  smooth.  Finally,  the  flows  down  the  landward-side  slope  in  the  small-scale  experiments 
had  little  if  any  air  entrainment.  This  is  certainly  not  the  case  for  similar  flows  at  full  scale.  How  this 
aeration  scale  effect  alters  the  parameters  estimated  by  the  empirical  equations  is  not  yet  known. 

ADDITIONAL  INFORMATION:  This  CHETN  is  a  product  of  the  Affordable  Levee  Strengthening 
and  New  Design  Work  Unit  being  conducted  at  the  U.S.  Army  Engineer  Research  and  Development 
Center,  Coastal  and  Hydraulics  Laboratory,  and  sponsored  by  the  Department  of  Homeland  Security 
Levee  Strengthening  and  Damage  Mitigation  Program.  Questions  about  this  technical  note  can  be 
addressed  to  Dr.  Steven  A.  Hughes  (Voice:  601-634-2026,  Fax:  601-634-3433,  email: 
Steven.A.Hughes@usace.army.mil).  Beneficial  reviews  were  provided  by  Norberto  C.  Nadal  and 
William  C.  Seabergh.  This  document  should  be  cited  as: 

Hughes,  S.  A.  2008.  Estimation  of  combined  wave  and  storm  surge  overtopping  at 
earthen  levees.  Coastal  and  Hydraulics  Engineering  Technical  Note  ERDC/CHL 
CHETN-III-78.  Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development 
Center,  http: //chi.  erdc.  us  ace.  army,  mil/chetn . 
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